Conservative chaotic map as a model of quantum many-body environment 
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We study the dynamics of the entanglement between two qubits coupled to a common chaotic 
environment, described by the quantum kicked rotator model. We show that the kicked rotator, 
which is a single-particle deterministic dynamical system, can reproduce the effects of a pure de- 
phasing many-body bath. Indeed, in the semiclassical limit the interaction with the kicked rotator 
can be described as a random phase-kick, so that decoherence is induced in the two-qubit system. 
We also show that our model can efSciently simulate non-Markovian environments. 

PACS numbers: 05.45.Mt,03.65.Yz,03.67.-a,05.45.Pq 



I. INTRODUCTION 

Real physical systems are never isolated and the cou- 
pling of the system to the environment leads to deco- 
herence. This process can be understood as the loss of 
quantum information, initially present in the state of the 
system, when non-classical correlations (entanglement) 
establish between the system and the environment. On 
the other hand, when tracing over the environmental de- 
grees of freedom, we expect that the entanglement be- 
tween internal degrees of freedom of the system is re- 
duced or even destroyed. Decoherence theory has a fun- 
damental interest, since it provides explanations of the 
emergence of classicality in a world governed by the laws 
of quantum mechanics [1]. Moreover, it is a threat to 
the actual implementation of any quantum computation 
and communication protocol [2, 3]. Indeed, decoherence 
invalidates the quantum superposition principle, which 
is at the heart of the potential power of any quantum 
algorithm. A deeper understanding of the decoherence 
phenomenon seems to be essential to develop quantum 
computation technologies. 

The environment is usually described as a many-body 
quantum system. The best-known model is the Caldeira- 
Leggett model [4] , in which the environment is a bosonic 
bath consisting of infinitely many harmonic oscillators 
at thermal equilibrium. More recently, first studies of 
the role played by chaotic dynamics in the decoherence 
process have been carried out [5, 6, 7, 8, 9]. These inves- 
tigations mainly focused on the differences between the 
processes of decoherence induced by chaotic and regular 
environments. 

In this paper, we address the following question: could 
the many-body environment be substituted, without 
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changing the effects on system's dynamics, by a closed 
deterministic system with a small number of degrees of 
freedom yet chaotic? In other words, can the complexity 
of the environment arise not from being many-body but 
from chaotic dynamics? In this paper, we give a positive 
answer to this question. We consider two qubits cou- 
pled to a single particle, fully deterministic, conservative 
chaotic "environment", described by the kicked rotator 
model. We show that, due to this system-environment 
interaction, the entropy of the system increases. At the 
same time, the entanglement between the two qubits de- 
cays, thus illustrating the loss of quantum coherence. We 
show that the evolution in time of the two-qubit entangle- 
ment is in good agreement with the evolution obtained 
in a pure dephasing stochastic model. Since this pure 
dephasing decoherence mechanism can be derived in the 
framework of the Caldeira-Leggett model [10], we have 
established a direct link between the effects of a many- 
body environment and of a chaotic single-particle envi- 
ronment. 

Finally, we show that also non-Markovian effects nat- 
urally appear in our model. This is interesting since such 
effects are relevant in several implementations of quan- 
tum information protocols [11], for instance with solid- 
state devices due to the coupling of the qubits to "slow" 
modes such as phonons. Moreover, non-Markovian ef- 
fects, apart from some exceptions [12], are usually very 
hard for analytic treatment. On the other hand, they 
can be included in models similar to the one discussed in 
the present paper with no additional computational cost 
with respect to the Markovian case. 

The structure of the paper is as follows: in Sec. II our 
model of a chaotic environment is introduced; in Sec. Ill 
both the entropy production and the entanglement decay 
are investigated, and the results compared with those ob- 
tained in the phase-kick model; in Sec. IV non-Markovian 
effects are analyzed; our conclusions arc drawn in Sec. V. 
The Bloch representation of the phase-kick map and the 
short time decay of the relevant bath (kicked rotator) 
correlation functions are discussed in Appendix A and 
Appendix B, respectively. 
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II. THE MODEL 

We consider two qubits coupled to a quantum kicked 
rotator. The overall system is governed by the Hamilto- 
nian 

H = + + ^(kr) ^ ^(int) ^ 

where H^^^ = oji (z = 1, 2) describes the free evolution 
of the two qubits, 

^('■■' = Y + fccos(^)^<5(r-jT) (2) 

3 

the quantum kicked rotator, and 

(int) ^ g (^(1) + ) cos(0) ^ (5(r - jT) (3) 

the interaction between the qubits and the kicked rota- 
tor; (Tq-* (a = x,y,z) denote the Pauli matrices of the 
«-th qubit. Both the cosine potential in ij(kr) ^^le in- 
teraction ijf'"*) are switched on and off instantaneously 
(kicks) at regular time intervals T. We consider the two 
qubits as an open quantum system and the kicked rotator 
as their common environment. Note that we chose non- 
interacting qubits as we want their entanglement to be 
affected exclusively by the coupling to the environment. 

The evolution from time tT (prior to the t-th kick) to 
time (t -I- f )T (prior to the {t + l)-th kick) of the kicked 
rotator in the classical limit is described by the Chirikov 
standard map: 

( n = n + ksinO, 

\e = e + Tn, (mod27r), ^' 

where (n, 0) are conjugated momentum-angle variables. 
Hereafter t = t/T denotes the discrete time, measured in 
number of kicks. By rcscaling n ^ p ^ Tn, the dynam- 
ics of Eqs. (4) is seen to depend only on the parameter 
K = kT. For A' = the motion is integrable; when 
K increases, a transition to chaos of the KAM type is 
observed. The last invariant KAM torus is broken for 
K K, 0.97. If /\ ~ 1 the phase space is mixed (simul- 
taneous presence of integrable and chaotic components). 
If K increases further, the stability islands progressively 
reduce their size; for K ^ \ they are not visible any 
more. In this paper, we always consider map (4) on the 
torus < < 27r, — tt < p < tt. In this case, the Chirikov 
standard map describes the stroboscopic dynamics of a 
conservative dynamical system with two degrees of free- 
dom which, in the fully chaotic regime K ^ \, relaxes, 
apart from quantum fluctuations, to the uniform distri- 
bution on the torus. 

The Hilbert space of the global system is given by 

n = n'^^^ ®n^'^'^ ®n^''''\ (5) 

where Ti^^^ and Ti'-^-' arc the two-dimensional Hilbert 
spaces associated to the two qubits, and TY^*-''') is the iV- 
dimensional space for the kicked rotator. 



The time evolution generated by Hamiltonian (1) in 
one kick is described by the operator 

exp [ - i(fc + e{ar^P + af^)) cos(^)] 

X cxp [ - iT^] exp ( - i (5i ai^^) exp {-182 ai^^) . 

(6) 

The effective Planck constant is fi-es = T = 2tt/N, where 
N is the number of quantum levels used to describe the 
kicked rotator; 5i = ujiT, 82 = i^2T; e is the coupling 
strength between the qubits and the environment. The 
classical limit hcff ^ is obtained by taking T — > and 
k — > 00, in such a way that K = kT is kept fixed. 

III. LOSS OF COHERENCE INDUCED BY A 
CHAOTIC ENVIRONMENT 

We are interested in the case in which the environment 
(the kicked rotator) is chaotic (we consider the kicked ro- 
tator with K ^ 1). The two qubits are initially prepared 
in a maximally entangled state, such that they are disen- 
tangled from the environment. Namely, we suppose that 
at t = the system is in the state 

1*0) = ® \%) , (7) 

where 10+) = ^ (|00) + |11)) is a Beh state (the par- 
ticular choice of the initial maximally entangled state is 
not crucial for our purposes), and jV'o) = X^n I") ^ 
generic state of the kicked rotator, with c„ random co- 
efficients such that J2n l^ril^ = and \n) eigenstates of 
the momentum operator. The evolution in time of the 
global system (kicked rotator plus qubits) is described 
by the unitary operator U defined in Eq. (6). Therefore, 
any initial pure state 1^*0) evolves into another pure state 
|5i(t)) = {7* l^'o). The reduced density matrix pi2(t) de- 
scribing the two qubits at time t is then obtained after 
tracing |5'(t)) (^(i)l o'^^r the kicked rotator's degrees of 
freedom. 

In the following we will focus our attention on the 
time evolution of the entanglement E12 between the two 
qubits and between them and the kicked rotator, mea- 
sured by the reduced von Neumann entropy S12. Clearly, 
for states like the one in Eq. (7), we have £^12(0) = 1, 
'S'i2(0) = 0. As the total system evolves, we expect that 
E12 decreases, while S'12 grows up, thus meaning that the 
two-qubit system is progressively losing coherence. 

The entanglement of formation E12 of a generic two- 
qubit state P12 can be evaluated following Ref. [13]. 
First of all we compute the concurrence, defined as 
C = max(Ai — A2 — A3 — A4,0), where the A^'s are the 
square roots of the eigenvalues of the matrix R = pi2Pi2, 
in decreasing order. Here pi2 is the spin flipped matrix of 
P12, and it is defined by pi2 = {(Jy®<Jy) P12 (cy^fy) (note 
that the complex conjugate is taken in the computational 
basis {|00), |01), |10), |11)}). Once the concurrence has 
been computed, the entanglement of formation is ob- 
tained as i? = /i((l -I- vT^~C^)/2), where h is the binary 
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entropy function: h{x 
with a; = (1 



)_= -x\0g2X-{l-x) l0g2(l -X), 

C^)/2. The other quantity stud- 



ied in our investigations is the entanglement between the 
two qubits and the kicked rotator, measured by the von 
Neumann entropy S = — Tr [pi2 log2 P12] of the reduced 
density matrix pi2. 

If the kicked rotator is in the chaotic regime and in 
the semiclassical region fieff ^ 1, it is possible to dras- 
tically simplify the description of the system in Eq. (1) 
by using the random phase-kick approximation, in the 
framework of the Kraus representation formalism. Since, 
to a first approximation, the phases between two consec- 
utive kicks in the chaotic regime can be considered as 
uncorrelated, the interaction with the environment can 
be simply modeled as a phase-kick rotating both qubits 
through the same random angle about the z-axis of the 
Bloch sphere. This rotation is described by the unitary 
matrix 



-ie cos 9 









-ie cos 9 









, (8) 



where the angle is drawn from a uniform random distri- 
bution in [0, 27r). The one-kick evolution of the reduced 
density matrix pi2 is then obtained after averaging over 



Pu^^J„d9R{e) 

xpi2 e»'^i'^i''e'''2<^<^)^t(6l). 



(9) 



In order to assess the validity of the random phase- 
kick approximation, we numerically investigate model (1) 
in the classically chaotic regime K ^ 1 and in the re- 
gion hcff ^ 1 in which the environment is a semiclassi- 
cal object. Under these conditions, we expect that the 
time evolution of the entanglement can be accurately 
predicted by the random phase model. Such expecta- 
tion is confirmed by the numerical data shown in Fig. 1. 
Even though differences between the two models remain 
at long times due to the finite number N of levels in the 
kicked rotator, such differences appear at later and later 
times when N 00 {h^s 0). The parameter K has 
been chosen much greater than one, so that the classi- 
cal phase space of the kicked rotator can be considered 
as completely chaotic. The actual value K sa 99.72676 
approximately corresponds to a zero of the Bessel func- 
tion J2{K). This is to completely wipe off memory ef- 
fects between consecutive and next-consecutive kicks (sec 
Sec. IV, Eqs. (13)-(14) for details). 

We point out that the random phase model can be 
derived from the Caldeira-Leggett model with a pure de- 
phasing coupling (X (ffi^^ -I- o-i^^) X]fc fffcft, with % cou- 
pling constant to the fc-th oscillator of the environment, 
whose coordinate operator is qk [10, 14]. This establishes 
a direct link between our chaotic single-particle environ- 
ment and a standard many-body environment. 

Numerical simulations show that, when the kicked ro- 
tator is chaotic, the pairwise entanglement E12 between 
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FIG. 1: Reduced von Neumann entropy 5i2 (main figure) and 
entanglement E12 (inset) as a function of time at if ~ 99.73, 
(5i = 10~^, 62 = y2(5i, e = 8 X 10"^. The thin curves corre- 
spond to different number of levels for the environment (the 
kicked rotator) (iV = 2^ 2^°, 2", 2" from bottom to 

top in the main figure and vice versa in the inset). The 
thick curves give the numerical results from the random phase 
model (9). 



the two qubits decays exponentially at short times, with 
a rate F proportional to e^: 



El2{t) 



(10) 



This is clearly seen from Figs. 2 and 3 [15]. In Fig. 2 
we plot the decay of entanglement E12 in time, for dif- 
ferent values of the oscillation frequencies 5i and 82 of 
the two qubits, at a fixed value of the coupling strength e 
with the environment. The short time decay of E12 is al- 
ways exponential and superimposed to oscillations whose 
frequency is determined by the internal dynamics of the 
two qubits [16]. Such oscillations can be clearly seen for 
5i = (52/\/2 = 10~^ (dashed curve), while their ampli- 
tude and period are too small to be seen on the scale of 
the figure for 5i = (52/V2 = 10^^ (solid curve). Finally, 
oscillations are absent at 5i = 82 /^/^ = 10~^ (dotted 
curve), as in this case their period is longer than the 
time scale for entanglement decay. Analogous remarks 
can be done about the oscillations in the entanglement 
production, measured by the reduced von Neumann en- 
tropy plotted in the inset of Fig. 2. 

An analytic estimate of the entanglement decay rate 
and of the entropy production can be derived from the 
random phase-kick model, in the limiting case (5i,(52 "C 
e ^ 1, that is, when the internal dynamics of the two 
qubits can be neglected on the time scale for entangle- 
ment decay. Starting from Eq. (9), it is possible to ex- 
plicitly write down the map pi2 — > pi2 in the Bloch rep- 
resentation, see Eqs. (A5) in App. A. Though in gen- 
eral these equations have to be solved numerically, in 
App. A we provide an explicit analytic solution valid for 
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FIG. 2: Entanglement E12 as a function of time for 
two qubits coupled to a chaotic kicked rotator {K ~ 
99.73) with A'' — 2^* levels and coupling strength 
e = 0.016. The various curves are for different os- 
cillation frequencies of the qubits; Si = 62 /V^ — 
10"^, (solid curve), 10~^ (dashed curve), 10~*(dotted curve). 
The thick dashed line gives the analytic estimate (11), valid 
in the regime 5i, S2 <^ e <^ 1. Inset: same as in the main fig- 
ure but for the Von Neumann entropy 512. The thick dashed 
curve is the estimate (12). 
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Eq. (10) for small times, with a « 4/ In 2. The analytic 
estimates (11) for E12 and (12) for S'12 are shown as thick 
dashed curves in Fig. 2. A very good agreement with 
numerical data can be clearly seen for ^2 ^ e <§; 1. 

In Fig. 3 the decay rate F of entanglement as a func- 
tion of e is shown; circles represent numerical data, while 
the dashed line is the analytic estimate of Eq. (11). A 
dependence F oc is found in both cases. In the inset we 
plot the characteristic time scale ts for the von Neumann 
entropy production. The time scale ts has been deter- 
mined from the condition £'12(^5) ~ C ~ 0.002 (the value 
chosen for C is not crucial). As before, circles stand for 
numerical data, while the dashed line is obtained from 
Eq. (12). From this figure we can conclude that the de- 
pendences of both F and ts on the system-environment 
coupling strength e are correctly reproduced by Eqs. (11) 
and (12), even though the oscillations due to the internal 
system's dynamics cannot be reproduced by these ana- 
lytic estimates. Finally, we point out that, as clear from 
Fig. 1, in the semiclassical chaotic regime {hcs <C 1 and 
^ 1) both F and ts can be reproduced with great 
accuracy by the random phase-kick model. 

We would like to stress that the results discussed in 
this section do not depend on the initial condition |-0o) 
in (7), provided that the kicked rotator is in the chaotic 
regime. On the other hand, we have found that both the 
entanglement decay and the entropy production strongly 
depend on lipo) in the integrable region K < 1. This 
implies that only in the chaotic regime a single particle 
can behave as a dephasing environment. 

Finally we notice that, contrary to other bath mod- 
els [14], the random phase-kick model does not generate 
entanglement between the two qubits. This has been nu- 
merically checked for model (1) both for initial pure sep- 
arable states and for separable mixtures. Moreover we 
also checked that, starting from a generic two qubit en- 
tangled state, the interaction with a chaotic memoryless 
environment cannot increase entanglement. Namely, we 
considered lO'' random initial conditions with i?i2(0) ^ 
and we found that, already after t ~ 100 kicks, entangle- 
ment has been always lowered: Ei2[t) < Ei2{0). 



FIG. 3: Decay rates F of the entanglement E12 (main figure) 
and time scales ts for the Von Neumann entropy (obtained 
from Si2{ts) = 0.002) (inset) as a function of the coupling 
strength eioi N = 2^^,K^ 99.73, Si = 82 /V^ = 10"^ Both 
numerical data (circles) and the analytic predictions (dashed 
curves) given by Eq. (11) (main figure) and Eq. (12) (inset) 
are shown. 

'^2 ^ e <C 1. In that case, for e^t ^ 1 we obtain: 

(11) 

Si2it) 0, ^{l-lni2e^t)), (12) 
consistently with the numerically obtained behavior of 



IV. NON-MARKOVIAN EFFECTS 

The model governed by Eq. (1) is very convenient for 
numerical investigations of decoherence. It can also be 
used to study physical situations which cannot be easily 
treated by means of analytic techniques using many-body 
environments like in the Caldeira-Leggett model. For in- 
stance, one could simulate, with the same computational 
cost, more complex couplings (e.g. when the direction of 
the coupling changes in time) or non-Markovian environ- 
ments. 

In this section, we show that memory effects naturally 
appear in our model, outside the range of validity of the 
random phase approximation. First of all, we point out 
that, given the coupling oc cos^, the bath corre- 
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FIG. 4: Oscillations of the entanglement decay rate F as 
a function of the classical chaos parameter K, at e — 5 x 
10-3, 5i = = 10-^ N = 2" (circles). Stars show 

the results obtained from a quantum trajectory approach 
(with average over 5000 trajectories) applied to the phase-kick 
model (16). Inset: entanglement decay rates as a function of 
for e = 5 X 10-^ 5i = S2 = 0, N = 2" (circles). The solid 
curve shows the analytic estimate (18). 



lation function relevant for the study of memory effects 
is {cos[9{t)]cos[9{t')]). Even in the chaotic regime, the 
phases are not completely uncorrelated: in the classical 
Chirikov standard map (4) correlations between cosines 
of the phases in two consecutive kicks are zero, but the 
same correlations between phases of two next-consecutive 
kicks do not vanish. As shown in Appendix B, we have 

(cos 6* cos 0) = 0, (13) 
(cos 61 cos I) = ^^^Y^, (14) 

where, given 6* = 9{t), 9 = 9{t+l) and 9 = 9{t+2), J2{K) 
is the Bcsscl function of the first kind of index 2 and K 
is the classical chaos parameter. Correlations between 
more distant kicks, (cos[0(i)] cos[6'(t')]), with t' — t > 2, 
though non vanishing, are very weak for K ^ 1, therefore 
hereafter we will neglect them. Correlations in the kicked 
rotator dynamics eventually result in a modification of 
the decay rate F of entanglement, as it is clearly shown 
in Fig. 4, where we plot the dependence of F as a function 
of K. 

In order to study analytically memory effects on the 
entanglement decay, we provide a simple generalization 
of the random phase-kick model (9), so that the corre- 
lations of Eqs. (13)-(14) are taken into account. We use 
the following conditional probability distribution for the 



angle 9 at time t + I, given the angle 9 at time t: 

f V^4^'-0+|(-l)*l +^^, (J2>0), 

p{9\9) = \ L 2 J 2n_ 

[ V^4^?-e-|]+i^^, (J2<0). 

(15) 

This distribution corresponds to having a 9 angle which 
is correlated with 9 (i.e. 9 ^ 9 ± tt/2) with probability 
Pc = ■\/ 1 J2 1 , and completely uncorrelated with probabil- 
ity pnc — 1 ~ Pc- For our purposes, the relevant prop- 
erties of the probability distribution p{9,9) = p{9)p{9\9) 
are: J d9 J d9p{9,9) = 1; J d9 J d9 p{9 , 9) cos 9 cos 9 = 0; 
/ d9 J d9p{9, 9) cos 9 cos 9 = J2/2. Therefore, Eq. (15) is 
a convenient distribution probability leading to a decay 
of the bath correlation function as in Eqs. (13)-(14) for 
the kicked rotator model. Of course, other probability 
distributions could equally well reproduce such decay. 

Given the probability distribution p{9\9), we can re- 
place (9) with the following two-kicks time evolution 
map: 

Pi2= d9 d9p{9,9)R{9)R{9)pi2R\9)R^{9), 
Jo Jo 

.(16) 

where pi2 and pi2 denote the two-qubit density matrix at 
times t and t+2 and p{9, 9) is the joint probability to have 
a rotation through an angle at a time t and through 
an angle 9 at time t + 1. Therefore, p{9,9) accounts 
for correlations between the angles at subsequent kicks. 
Clearly, if angles are completely uncorrelated we have 
p{9,9) = l/47r^, thus recovering the random phase-kick 
model (9). Note that, for the sake of simplicity, map (16) 
has been written for 61 = 62 = (the generalization of 
this map to (5i , is straightforward) . 

The phase-kick model (16) can be simulated by us- 
ing the quantum trajectories approach [17, 18, 19]. This 
method is very convenient in the study of dissipative sys- 
tems: instead of solving a master equation, one stochasti- 
cally evolves a state vector, and then averages over many 
runs. At the end, we get the same probabilities as the 
ones directly obtained through the density matrix. In our 
case, the effect of the kicked rotator on the two-qubits 
wave function is simply that of a rotation through an 
angle whose value is drawn according to the probabil- 
ity distribution (15). Numerical data obtained with the 
quantum trajectories method are plotted in Fig. 4 (stars); 
notice that, a.t K 1, they are in good agreement with 
data from simulation of the Hamiltonian model (1) (cir- 
cles). 

It is also possible to give an analytic estimate of the 
decay rate T{K) in the limit in which the free evolu- 
tion of the two qubits can be neglected (i.e., we take 
Si,S2 <C e ^ 1). As for the random phase model (9), in 
this limit the effect of map (16) is pure dephasing. At 
every map step the density matrix pi2 is of the form of 
Eq. (A8). The coherences Dt can be computed by iter- 
ating Eq. (16). We obtain what follows for the first time 
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steps: 

= l-4e2, 

D3 = l-[l2 + 8J2{K)y, ^^'> 
Di = l- [l6+16J2{K)]e\ 

Assuming an exponential decay of entanglement, 
Ei2{t) ~ e~'"*, we can evaluate T starting from D3 and 
D4: r « ln(£;3/£:4). Thus, from Eqs. (17) we obtain the 
following analytic estimate: 

r« ^(1 + 2,72 W). (18) 

In the inset of Fig. 4 we compute the entanglement 
decay rates as a function of K in the case 5i — 82 = 0. 
This figure clearly shows that the rates obtained from nu- 
merical data for the chaotic environment model (circles) 
are in good agreement, when ilT 3> 1, with the analytic 
estimate (18) (solid curve). 

We point out that the oscillations with K of the en- 
tanglement decay rate T{K) are ruled by the Bessel func- 
tion J2{K), in the same way as for the well-known K- 
oscillations of the classical diffusion coefficient D{K) [20]. 
Therefore, the entanglement decay rate is strictly related 
to a purely classical quantity. The ultimate reason of 
this relation is rooted in the fact that the bath corre- 
lation function relevant for the study of memory effects, 
{cos[d{t)] cos[9{t')]), also determines the deviations of the 
diffusion coefficient D{K) from the random phase ap- 
proximation value K"^ /2. 



V. CONCLUSIONS 

We have shown that a single particle in the chaotic 
regime can be equivalent to a bath with infinitely many 
degrees of freedom as a dephasing environment. This pro- 
vides a clear example of the relevance of chaos to dccoher- 
ence. Furthermore, also memory effects can be included 
in a single particle chaotic environment. This observa- 
tion paves the way to a convenient numerical simulation 
of non-Markovian dynamics. Indeed, both exponential 
and power-law decays of the bath (single particle) corre- 
lation functions can be obtained in chaotic maps [21, 22]. 
Therefore, these maps could be used to efficiently simu- 
late important quantum noise models such as random 
telegraphic or 1// noise. 
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APPENDIX A: THE RANDOM PHASE MODEL 

Here we provide an explicit set of equations describing 
map (9) in the Bloch representation. In such representa- 
tion, a generic two-qubit mixed state can be written as 

^ z 3 3 

Pi2 = -4.®l+^a^{a,®l)+^l3j{l®aj)+^ lij{oi®aj) 

i=l j=l = l 

(Al) 

We insert this expansion into Eq. (9) and evaluate the 
commutators between pi2 and the terms that multiply 
P12 on the right in (9) (e^'^i'^-", e**'^"^-", and R'^{9)). For 
this purpose, we use the standard commutation rules for 
the Pauli matrices: 

[(Ji,(Tj] = 2ieijkak , (A2) 

where e^-fc is the Levi-Civita tensor, which is equal to 
-|-1 if ijk is an even permutation of 123, to —1 for odd 
permutations and otherwise. We also use the identity 

^^sa^ = (cos S)i + i (sin S) (A3) 

and expand the exponentials in W {9) to the second order 
in e. The average over 9 is then performed using 

(cos2"+i6')„ = 0, 

^ (A4) 
(cos^e), = l, (cos^e),^|. 

In this way we arrive at the following one-kick map, 
valid in the approximation e, (5i, (^2 <C 1: 

ctx = ax{l - e^), 
Uy = ay{l - e^) - 2Siaz, 
Q-z = az + 2Siay, 



py=Py{l-e^)-2S2f^z, 

Pz=f3z + 2S2Py, 





= lxx{^ 


-2e2) 


+ 2e^7TO, 


%y 


= 7to(1 


-2e2) 


+ 26^7:^:1, - 2(5i7^y - 252lyz, 


Izz 


= Izz + 


25i7i/z 


+ 2(527^y7 


Ixy 


= lxy{'^ 


-2e2) 


- 2e'^-fyrc ~ 252lxz, 


"fyx 


= lyx{l 


-2e2) 


- 26^7^^ - 282! zx, 


Ixz 


= lxz{l 


-2e2) 


+ 2^27x1;, 


Izx 


= lzx{l 


-2e2) 


+ 25i7ya;, 


lyz 


= 7y^(l 


-2e2) 


- 2(5i722 -f 252lyy, 


. Izy 


= lzy{^ 


-2e2) 


H- 2(5i7.y.y - 252'-izz- 



(A5) 

Note that these equations can be straightforwardly gen- 
eralized to the continuous time limit. 
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We now focus on the case in which the two qubits are 
initially in the Bell state | </>"''), so that the only non zero 
coordinates in the Bloch representation of puit = 0) = 
{(f>+\ are j^x = Izz = -lyy = 1/4. The system in 
Eq. (A5) then reduces to a set of five coupled equations 
for the 7's: 



lyy = 7y;/(l - 2e2) + 26^7^^ - 2(5i7^y ■ 

lyz = lyz{^ - 26^) - 2(5i722 + 2(527yy, 
Izy = 72y(l - 26^) + 25i-fyy - 2(527zz- 



2^271/2 1 



(A6) 

Their solution after t map steps is given, in the limiting 
case 61,82 < e < 1, by 



lxx{t) ~ -lyyit) 



e\ lzz{t) 



(A7) 



Note that here we also require e^t <C 1. The correspond- 
ing density matrix is 



with 



1 L»t 





L>t 1 



(A8) 



(A9) 



The calculation of the von Neumann entropy Si2{t) and 
of the entanglement Ei2{t) for state (A8) is straightfor- 
ward and leads to Eqs. (11)-(12). 



APPENDIX B: ANGULAR CORRELATIONS IN 
THE CHIRIKOV STANDARD MAP 



where K = kT. In this appendix we provide a simple 
proof of Eqs. (13)-(14), namely we evaluate correlations 
between the cosines of the phases of two consecutive and 
next-consecutive kicks. Averages are performed over all 
the input phase space: 



{f{e,p)) EE hm 

P—>^OG 




(B2) 



By using Eq. (Bl), we obtain the following equalities for 
correlations between two consecutive kicks: 

(cos 9 cos 9') ~ (cos 9 cos {9 + p + K sin 9)) 

= (cosp) (cos 6* cos (6* -I- A' sin 6*)) (B3) 

- (sinp) (cos e sin [9 + K sin 6*)) = 0, 
where we have used some trivial trigonometric identities, 
and (cosx) = (sinx) = 0. 

Correlations between two next-consecutive kicks are 
evaluated by considering values at times t + I (denoted 
by an overbar) and t — 1 (denoted by an underbar). The 
map in Eq. (Bl) can be straightforwardly inverted to give 
backward evolution: 



p = p — K sm{9 — p), 
1^9-p, (mod 27r). 



(B4) 



Therefore, from Eqs. (Bl) and (B4), after some simple 
algebra we obtain: 



(cos 6' cos ^) = (cos {9 + p + K sm9) cos {9 — p)) 

= 1 (cos (20 + A- sine)) = \j-2{K) - \j2{K), ^^^^ 



where 



After rescaling the angular momentum n 
the Chirikov standard map (4) becomes 



p + K sin 9, 

9+p, (mod 27r), 



p = Tn, 



(Bl) 



1 r 

Jn (K) = — / cos (a: sin 6* - n9) d9 (B6) 
27r Jo 



is the Bessel function of the first kind of index 2. 
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